Comparative analysis of tunneling magnetoresistance in low- Tc Nb/AlAlOx/Nb and 
high-Tc Bi2 yPbySr2CaCu208+5 intrinsic Josephson junctions 

V. M. KrasnovQ H. Motzkau, T. Golod, A. Rydh, and S.O. Katterwe 
Department of Physics, Stockholm University, AlbaNova University Center, SE-10691 Stockholm, Sweden 

A.B. Kulakov 

Institute of Solid State Physics, Russian Academy of Sciences, 142432 Chernogolovka, Russia 

(Dated: January 20, 2013) 

We perform a detailed comparison of magnetotunneling in conventional low- Tc Nb/AlAlOx/Nb 
junctions with that in slightly overdoped Bi2-j/PbySr2CaCu208+6 [Bi(Pb)-2212] intrinsic Josephson 
junctions and with microscopic calculations. It is found that both types of junctions behave in a 
qualitatively similar way. Both magnetic field and temperature suppress superconductivity in the 
state-conserving manner. This leads to the characteristic sign-change of tunneling magnetoresis- 
tance from the negative at the sub-gap to the positive at the sum-gap bias. We derived theoretically 
and verified experimentally scaling laws of magnetotunneling characteristics and employ them for 
accurate extraction of the upper critical field Hc2- For Nb an extended region of surface super- 
conductivity at Hc2 < H < Hc3 is observed. The parameters of Bi(Pb)-2212 were obtained from 
self-consistent analysis of magnetotunneling data at different levels of bias, dissipation powers and 
for different mesa sizes, which precludes the influence of self-heating. It is found that Hc2{0) for 
Bi(Pb)-2212 is ~ 70 T and decreases significantly at T ^ Tc. The amplitude of sub-gap magnetore- 
sistance is suppressed exponentially at T > Tc/2, but remains negative, although very small, above 
Tc- This may indicate existence of an extended fiuctuation region, which, however, does not destroy 
the general second-order type of the phase transition at Tc. 

PACS numbers: 74.55.-l-v 74.25.0p 74.72.Gh 74.50.-|-r 



I. INTRODUCTION 

Magnetoresistance (MR) is one of the basic tools for 
analysis of the electronic structure of metals. MR data 
contains important information about bulk electronic 
structure of cuprate high temperature superconductors. 
This was demonstrated by recent MR studies providing 
compelling evidence for reconstruction of the Fermi sur- 
face in the pseudogap state of underdoped cuprates as a 
result of density wave ordering [J Q • Several experimen- 
tal techniques revealed the existence of two distinct en- 
ergy scales in cuprates: the superconducting gap A and 
the normal state pseudogap (PG) (for review, see e.g. 
Refs. [3,-^6,]) with different behavior with respect to tem- 
perature d, ItI-ITsI , doping [l|, [13, [11] and magnetic field 
[l6l-(l8l| . However, there is no consensus on whether they 
are competing [ij, [l^, cooperating Q, or representing 
two manifestations of the same phenomenon [1,H, [1^ [2l| . 

Analysis of magnetic field effects is particularly use- 
ful for scrutinizing the superconducting origin of the 
gaps. Non-superconducting (e.g., structural, antiferro- 
magnetic, charge, spin, or d-density wave) orders are 
typically insensitive to achievable fields [l^ . In this case, 
magnetic field may selectively suppress the superconduct- 
ing gap. However, discussion of magnetic field effects in 
cuprates remains controversial. Confiicting reports ex- 
ist even on such a basic parameter as the upper critical 
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field Hc2 ■ It was reported that Hc2 behaves in a conven- 
tional manner, i.e., vanishes at Tc and scales with Tc as a 
function of doping [22 - 26]. But it was also reported that 
Hc2 is T-independent and persists well above Tc plj and 
increases with underdoping despite reduction of Tc [13] . 

There are many obstacles for deciphering MR data of 
cuprates such as: fuzzy superconducting transition due 
to persistence of the PG; d-wave symmetry of A; ill de- 
fined quasiparticles (QP) and strong angular dependence 
of QP scattering rates [28|; very high anisotropy, which 
requires accurate control of transport current direction 
and sample geometry; existence of an extended fluctua- 
tion region above Tc [29| ; possible doping inhomogeneity 
[^; and extremely large Hc2 ~ 100 T. One of the im- 
portant open question is whether so large magnetic fields 
simply suppress superconductivity, or simultaneously in- 
duce a competing order, as may be suggested by obser- 
vation of charge density and antiferromagnetic spin 
order [s^] in vortex cores. 

Typically MR involves only QPs at the Fermi level, 
averaged over the Brillouin zone, and thus does not pro- 
vide spectroscopic information about the QP density of 
states (DOS) away from the Fermi surface. A rare ex- 
ception is the c-axis transport in extremely anisotropic 
layered cuprates. Single crystals of Bi, Tl [3l|, [sl] 
and Hg [Sj] based cuprates represent natural stacks of 
atomic scale intrinsic tunnel junctions. The intrinsic tun- 
nelling spectroscopy (ITS) provides a unique opportu- 
nity to probe directly bulk electronic spectra of cuprates 
[IM [ll, [11 [3. Tunneling MR is potentially a very 
powerful tool for analysis of superconducting features in 
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electronic spectra. This was demonstrated in previous 
studies for cuprates H, [H, [35l - l37| as well as for conven- 
tional low- Tc [H-Iijand non-cuprate [Hj high Tc su- 
perconductors. Application of magnetic field leads to 
appearance of a spatially inhomogeneous mixed state. 
This makes analysis of magnetotunneling data non-trivial 
[4^ and even counterintuitive [3] ■ Therefore, a clear 
understanding of how the magnetotunneling characteris- 
tics of Josephson junctions should behave is needed for 
accurate data analysis. 

In this work we perform detailed comparison of magne- 
totunneling in low-Tc Nb/AlAlOx/Nb, and slightly over- 
doped Bia-yPbj^SrsCaCuaOg+a [Bi(Pb)-2212] intrinsic 
Josephson junctions with theoretical calculations. Small 
sizes of our Bi(Pb)-2212 mesas, which are one-two or- 
ders of mag nitude smaller than in previous similar studies 
(igI ITsI . [36| . in combination with the ability to extract in- 
formation at the sub-gap bias with low dissipation power, 
lead to effective obviation of self- heating IH, l46| . 
Both low- and high- Tc junctions show qualitatively sim- 
ilar behavior. Magnetic field and temperature suppress 
superconductivity in the state-conserving manner: en- 
hancement of the sub-gap conductance due to suppres- 
sion of A(T, H) is exactly compensated by reduction of 
the sum-gap conductance peak. As a result, the MR 
changes sign from the negative at the sub-gap voltages 
to the positive at the sum-gap eV > 2A. This allows 
us to trace closing of A at T ^ Tc with unprecedented 
clarity. We derive simple scaling laws for magnetotunnel- 
ing and employ them for unambiguous extraction of su- 
perconducting parameters such as A{T,H) and Hc2{T). 
The extracted Hc2{T) decreases significantly upon ap- 
proaching Tc. Our data indicates that superconductivity 
in slightly overdoped cuprates appears in a conventional 
manner by means of the second order phase transition. 

The paper is organized as follows. In sec. II we de- 
scribe the theoretical formalism, used for microscopic cal- 
culation of current- voltage (/ — V) characteristics in the 
mixed state. Sec. Ill describes the experimental setup 
and studied junctions. Methods used for obviation of 
self-heating are discussed. In sec. IV we present main 
results. It is shown that both Nb and Bi-2212 junctions 
behave in a qualitatively similar way as a function of T 
and H. In sec. V we analyze scaling laws of different 
tunneling parameters as a function of H/Hc2 and apply 
them for extraction of Hc2 in Nb-j unctions. This also 
provides a clear evidence for persistence of the surface 
superconductivity in Nb. In sec. VI we apply the same 
scaling laws for extraction of Hc2{T) for Bi(Pb)-2212 and 
analyze the remaining MR above Tc. Finally, in sec. VII 
we summarize our conclusions. 



II. THEORY 

To calculate spatially non-uniform distribution of the 
gap and the DOS in the mixed state, we use the circular 
cell approximation [47l - l49| . We approximate the hexago- 



nal unit cell of the Abrikosov vortex lattice by a circular 
cell with a radius corresponding to one flux quantum $0 
within the cell, see the sketch in Fig. [T] (a). 
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and assume a cylindrically symmetric vector potential 

Q{r) = l/r-r/p'. (2) 

Spatial distribution of the order parameter A(r) at dif- 
ferent temperatures and magnetic fields are calculated by 
solving microscopic Usadel equations: 
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g2(r) sin(e„) cos(e„) + A cos(e„) = 0, 
together with the self-consistency equation, 

A ln(T/Tc) + 2ttT ^ ( A/w„ - sin(e„)) = 0. (4) 

n 

Here ^5 = {D /2tiTc)'^/'^ , where D is the QP diffusion 
coefficient, lli„ = 7TT(2n+ l),n = 0,1,2,... are Matsub- 
ara frequencies, cos(0„) and sin(0„) are the normal and 
the anomalous Green function components, and primes 
denote spatial derivation d/dr. Those equations are sub- 
ject to boundary conditions at the center of the vortex 
r = 0: A(0) = e„(0) = 0, and at the edge of the cir- 
cular cell r = p: /S.{p)' = Qn[p)' = 0. This system of 
nonlinear equations was solved numerically, up to a cut- 
off frequency n ~ 50 T^/T ^ using an iterative procedure. 
More details about the formalism and the numerical pro- 
cedure can be found in Refs. [isl - fsoj . 

Fig. [1] (b) shows an example of calculated spatial 
distribution of A(r) along the line connecting two vor- 
tices for the case of Nb (Tc = 8.8iC) at T = 4,.7K and 
H/Hc2 = 0.5. The A(r) is normalized by the equilibrium 
value of the gap Ao(T) at H = 0. 

To calculate I — V characteristics in the mixed state, 
we also need to calculate the spatial variation of the QP 
DOS N{r, E) as a function of the QP energy E. This is 
done by analytic continuation of discrete 0„(a;„) to the 
continuous energy axis via substitution, a;„ = —iE, in 
Eq. ([3]), where i is the imaginary unit: 
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iEsm{Q)~ (5) 

7 J 

f{r) sin(e) cos(e) + A(r) cos(e) = 0. 



The spatially non-uniform DOS is then obtained as 
N{r,E) = 5R(cos[0(r, £;)]). 



(6) 



Lines in Fig. [T] (c) show N{r,E), normalized to the 
DOS in the normal state, at three points of the circular 
cell: A - at the center of the vortex, B - at the inter- 
mediate point where A(r) recovers to half of its maxi- 
mum value, see Fig. [1] (b), and C - at the edge of the 
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FIG. 1. (Color online), (a) A sketch of the hexagonal Abrikosov vortex lattice and the circular approximation of a unit cell, 
(b) Calculated spatial distribution of the order parameter between two vortices, (c) Local density of states at three points A, 
B, C, indicated in (a) and (b). Calculations are made for H = 0.5 Hc2, T — 4.7 K and Tc = 8.8 K. 



cell. In the center of the vortex the superconductivity 
is completely suppressed, A = 0, and N{E) = 1, as in 
the normal state. Away from the vortex the order pa- 
rameter is partly restored, but QP spectra are gapless 
N{E < A) 7^ di^iO]. Due to spatial inhomogeneity, 
the maximum in DOS is significantly smeared, compared 
to the BCS singularity, and the energy of the maximum 
is no longer equal to A(r). 

Tunneling I — V characteristics are calculated by inte- 
gration over the circular unit cell: 



1 f 2r 

Rn Jo P J-oo 

N{r, E)N{r, E -t- eV)f{E) [1 - f{E + eV)] . 



(7) 



Here i?„ is the tunneling (normal) resistance and / is the 
Fermi distribution function. 



A. Validity of the model 

The described formalism is valid for arbitrary T and H 
for dirty type-II superconductors with s-wave symmetry 
of the order parameter, which is appropriate for sput- 
tered Nb films. Calculations presented below are made 
for BCS parameters with T^. = 8.8 K, typical for Nb. 
This allows direct comparison with experimental data for 
Nb/AlAlOx/Nb junctions. 

In Eq. ([7]) we disregard possible misalignment of vor- 
tices in the two electrodes. For intrinsic Josephson junc- 
tions, due to very high anisotropy and weak interlayer 
coupling, such misalignment may be significant at low 
magnetic fields ^] , but could be neglected for high fields 
used in this work. 

Certain deviations can be expected for Bi-2212 due 
to the d-wave symmetry of the order parameter, which 
reduces the sum - gap singularity and affects the vor- 
tex structure |5ll-l53j. However, theoretical calculations 



demonstrated [42|, |5l| that the scaling of DOS charac- 
teristics as a function of H/ Hc2 , which will be discussed 
below, is valid also for d-wave superconductors (see e.g. 
Figs. 5 from Ref. [5l|). Therefore, the procedure of ex- 
traction of Hc2 from such a scaling should be also valid 
for intrinsic Josephson junctions. 

Coexistence of competing order parameters, associated 
with the pseudogap, could also affect the tunneling DOS 
in cuprates [13, H^- An influence of the PC on the intrin- 
sic tunneling MR was reported for underdopcd Bi-2212 
[S E, [13, Si- Namely, with underdoping the sub- 
gap MR is rapidly reduced [27[ and suppression of su- 
perconductivity by magnetic field occurs in a seemingly 
non-state-conserving manner [36j . This may indicate a 
gradual recovery of the competing PG order upon sup- 
pression of superconductivity, as observed in the vortex 
cores [1, [S^l- To avoid possible complications, we restrict 
our analysis to slightly overdoped Bi(Pb)-2212, for which 
there is no significant distortion by the PG [T3ll56|. 



III. EXPERIMENTAL 

We study tunneling magnetorestistance in standard 
low-Tc Nb/AlAlOx/Nb junctions and in small slightly 
overdoped Bi(Pb)-2212 mesa structures containing few 
atomic scale intrinsic Josephson junctions. Measure- 
ments were performed in a gas-fiow ^He cryostat in a tem- 
perature range down to 1.6 K and magnetic field H up to 
17 T. Samples were mounted on a rotatable sample holder 
with the alignment accuracy better than 0.02°. Details 
of the measurement setup can be found in Refs. [s^. [57|. 



A. Nb/AlAlOx/Nb junctions 

Nb/AlAlOx/Nb junctions were made by the standard 
HYPRES trilayer technology ^8| with a critical current 
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density of 1000 A/cm^. A detailed description of junction 
parameters can be found in Ref . [59|] . Junctions consist of 
two sputtered Nb thin films with thicknesses 150 and 50 
nm for base and counter electrodes, respectively. Due to 
different thicknesses, electrodes have slightly different Tc- 
The critical temperature of the junction is ~ 8.8 K. The 
junction barrier is formed by deposition of a thin Al layer 
with the thickness ^ 10 nm on top of the base electrode, 
followed by a subsequent oxidation to form the AlOx tun- 
nel barrier. During oxidation, only the surface layer of 
Al is oxidized, leaving the rest of Al intact. This results 
in a proximity effect between the bottom Nb layer and 
Al [50| ■ Tunneling occurs between the proximity-induced 
superconducting layer of Al and the top Nb layer. A de- 
tailed analysis of the proximity effect in Nb/AlAlOx/Nb 
junctions can be found in Ref. [soj . 

Several junctions with different sizes on the same chip 
were studied and showed similar results. As for the case 
of Bi-2212 mesas [l3|j with increasing junction area the 
sum- gap kink m. I—V becomes excessively sharp and may 
even develop a back-bending as a result of progressive 
self- heating Self-heatin g is effectively obviated by 

miniaturization of junctions |l3l l45l l46j . Therefore, in 
what follows we show data only for the smallest junction 
with sizes ^ 2.5 x 2.5 /xm^, which is least affected by 
self-heating. 



B. Bi-2212 intrinsic Josephson junctions 



over-doping [l^ . 

The c-axis transport in Bi-2212 is non-metallic due to 
interlayer tunneling mechanism of transport, in combi- 
nation with the so-called c-axis pseudogap. The latter 
is much morepronounced than the PG in the a6-plane 
transport [1, U HI] and exists in a broader temperature 
0, [iJI and doping ranges @ . In Ref. [3J] it was shown 
that in the normal state T > Tc, c-axis intrinsic tunneling 
characteristics exhibit a trivial thermal-activation behav- 
ior, described by just one constant - the effective barrier 
height. This puts a question whether the c-axis PG repre- 
sents the real two-particle gap in the DOS or is just the 
single QP tunneling matrix phenomenon. In Ref. (Tsj 
more subtle features in intrinsic tunneling characteristics 
were found beyond the thermal-activation background. 
Those appear in the same temperature region as in the 
a&-plane transport and were attributed to the genuine 
two-particle pseudogap in the QP DOS. A similar con- 
clusion about the existence of two distinct pseudogap-like 
phenomena above was also reached in optical infrared 
ellipsometry studies [1,16^1 • The PG complicates analysis 
of tunneling magnetoresistance because it makes the su- 
perconducting transition fuzzy. Since the main purpose 
of this work is to establish how to extract unambiguous 
information out of magnetotunneling data, here we will 
mostly concentrate on analysis of overdoped Bi(Pb)-2212 
crystals, which according to previous studies are less af- 
fected by the PG and are most close to the conventional 
BCS-type superconductivity [scj . 



We study small, micrometer-size mesa structures con- 
taining few atomic scale intrinsic Josephson junctions. 
The mesas are fabricated on top of Bi-2212 single crys- 
tals using micro/nano-fabrication techniques. Details of 
the cryst al g rowth and sample fabrication can be found 
in Refs. j6lj and [13], respectively. Several contacts on 
top of the crystals allow us to perform three or quasi-four 
probe measurement of the pure c-axis transport [63| . De- 
tails of measurements and mesa characterization can be 
found elsewhere [ll,[il|3l]. 

We present data for two batches of crystals: lead- 
doped, slightly overdoped Bi2-j^Pbj^Sr2CaCu208-(-x with 
Tc ~ 89 — 93 K and yttrium-doped, slightly underdoped 
Bi2Sr2Cai_:rY:,Cu208+:, (Bi(Y)-2212) with ~ 92 K. 
It should be said that the c-axis phase coherence in 
small mesas is not a good measure of Tc- The asso- 
ciated Josephson coupling energy density is small, and 
the total energy is decreasing proportional to the mesa 
area. Therefore, in small mesas the phase coherence 
and the c-axis critical current are suppressed by ther- 
mal fluctuations [12, at temperatures significantly 
lower than that for the in-plane transport. A more de- 
tailed discussion on determination of Tc can be found in 
Ref. [l^. Both types of crystals have similar optimal 
Tc ~ 96 K. The most noticeable difference between them 
is in the c-axis critical current density Jc(4.2K) ~ IC* 
and lO^A/cm^ for lead- and yttrium-doped mesas, re- 
spectively. This is due to a rapid increase of Jc with 



C. Obviation of self-heating in Bi-2212 mesas 

Mesas with different sizes were made on the same single 
crystal. As in the case of Nb/AlAlOx/Nb junctions |60|, 
we observe that I — V characteristics of larger mesas are 
more distorted by self- heating at large bias [13[. To ob- 
viate self-heating, we perform additional miniaturization 
of mesa structures [13i, i45i , ,46 ] down to sub-micrometer 



sizes using Focused Ion Beam trimming 62 1. 

The dissipation powers at the sum-gap knee va I — V 
(the peak in dl/dV) for the two small mesas shown in Fig. 
[1(c) and (d) are P ~ 0.13 and 0.19 mW for Bi(Pb)-2212 
and Bi(Y)-2212 mesas, respectively, at the lowest tem- 
peratures. According to previous reports [H,!!^, typical 
thermal resistances of our micrometer-size mesas lie in 
the range ~ 70 — 30 K/mW at low T and decrease to ^ 10 
K/mW at Tc- Thus, self- heating at the sum-gap peak for 
those small mesas is manageable and the effective mesa 
temperature at the peak remains well below Tc- For the 
Bi(Y)-2212 mesa this was unambiguously proven by anal- 
ysis of the size-dependence of intrinsic tunneling spectra 
for mesas with different sizes on the same single crystal 
[H, HI] . To completely exclude possible artifacts of self- 
heating from the data analysis, in what follows we define 
superconducting parameters from scaling laws valid for 
tunneling magnetoresistance at different bias levels, in- 
cluding low, sub-gap bias with negligible heating. 
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FIG. 2. (Color online). Temperature dependence of tunneling I ~ V characteristics at zero magnetic field for (a) calculations 
for a BCS superconductor (Nb, Tc = 8.8 K), (b) for the Nb/AlAlOx/Nb junction, (c) for slightly overdoped Bi(Pb)-2212 mesa 
and (d) for slightly underdoped Bi(Y)-2212 mesa, (e) and (f) Measured T-dependence of the sum-gap voltage. Solid lines 
represent BCS T-dependence. (g) and (h) Comparison of the correlation between the sub-gap resistance (solid symbols) and 
the sum-gap conductance (open symbols) for the same Nb/AlAlOx/Nb junction and Bi-2212 mesas, respectively. 



IV. RESULTS AND DISCUSSION 

A. Temperature dependence of tunneling 
characteristics at zero magnetic field 

Figure [5] shows T-dependencies oi I — Vs a.t H — 
for (a) theoretical calculations for Nb parameters, (b) 
the Nb/AlAlOx/Nb junction, (c) a slightly overdoped 
Bi(Pb)-2212 mesa with the in-plane area 0.9 x 1.3 firn^ 
and Tc ~ 91 K and (d) a slightly underdoped Bi(Y)-2212 
mesa with the area 1.8 x 2.0 ^m^ and Tc ~ 92 K. Both 
mesas contain = 9 intrinsic Josephson junctions, es- 
timated by counting QP-branches in I — V [13]. It is 
seen that / — T^s of Bi-2212 mesas are closely resem- 
bling those for conventional superconductor-insulator- 
superconductor tunnel junctions. The pronounced sum- 
gap kink is clearly seen at low temperatures, followed by 
an almost T-independent tunneling resistance at higher 
bias. The sum-gap kink moves to lower voltage and van- 
ishes in amplitude upon approaching the Tc [3, [Tsj . 

Experimental characteristics of the Nb/AlAlOx/Nb 
junction can be explicitly compared with the correspond- 
ing theoretical simulations. It is seen that there is a good 
agreement, however, the proximity effect between Al and 
Nb leads to some smearing out of the sum-gap kink and 
enhancement of the sub-gap conductivity. The most clear 
proximity induced peculiarity in / — of Nb/AlAlOx/Nb 
is the pronounced dip in dl/dV above the sum-gap peak 



(see Fig. IHb). The dip is caused by the double max- 
ima structure of the DOS in Al ^50^] : the lower maximum 
corresponds to the proximity induced energy gap in Al, 
A^;, the upper - to the inherited gap from Nb, Ajvb. 
In this case the sum-gap peak in conductance occurs at 
eVsg — Aai + A^vb and the dip at eV ~ 2Ajvb. Thus 
determined gaps are: Aai = 1.22 meV and A^i, = 1.49 
meV, consistent with previous reports |5fl]. 

Due to the d-wave symmetry of the order parameter in 
cuprates, the sum-gap kinks in Bi-2212 intrinsic tunnel- 
ing characteristics are significantly more smeared than in 
s-wave low- Tc junctions. From comparison of Figs. [2] (c) 
and (d) it is also clear that the sum-gap kink is further 
loosing its sharpness with underdoping [14.] . 

Figure [2] (e) and (f ) represent T-dependencies of the 
sum-gap voltage (peak in dl/dV) at H = for the 
Nb/AlAlOx/Nb junction and the Bi-2212 mesas (per in- 
trinsic Josephson junction), respectively. They were ob- 
tained from I — V shown in Figs. [2] (b-d). The bulk gap 
of Bi(Y)-2212 follows very accurately the standard BCS 
dependence. More details on the T— dependence of inter- 
layer tunneling characteristics of our mesas can be found 
in Refs. m,^]. 

Both Nb and Bi-2212 demonstrate signatures of strong 
coupling superconductivity. For Nb the ratio 2Ajv(,/Tc ~ 
3.9, larger than the weak coupling BCS value of 3.5 for 
s-wave superconductors. For the slightly underdoped 
Bi(Y)-2212, 2Ao/Tc ~ 8.2. In the slightly overdoped 
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FIG. 3. (Color online). Magnetic field evolution of / — V characteristics at low T for (a) BCS calculations at T ~ 2 K and 
Tc = 8.8 K, (b) the Nb/AlAlOx/Nb junction at T ~ 2 K. Semi-logarithmic plots of (c) theoretical 7 - Fs at T ~ 0.5 K, and 
(d) J — Vs of a Bi(Pb)-2212 mesa in the sub-gap region at T ~ 2 K for different c-axis field components. Note a remarkable 
almost parallel translation of / — V^s in the semi-logarithmic scale in (c) and (d). 



Bi(Pb)-2212, the ratio decreases to 6.1. This is consis- 
tent with previous break junction studies, which indi- 
cated that in overdoped Bi-2212 2Ao/Tc is approaching 
the weak coupling BCS value for d-wave superconductors 
of 4.28 [H. 

As seen from Figs. [2] (b-d) in all cases the increase of 
T leads to simultaneous increase of the sub-gap current 
and decrease of the sharpness of the sum-gap kink. To 
demonstrate this reciprocal correlation, in Figs. [l](g) and 
(h) we plot the T— dependence of the sub-gap resistance, 
dV/dI{V = A/e), (solid symbols) and the sum-gap con- 
ductance peak, dI/dV{V = 2A/e), (open symbols) for 
the same junctions. For Nb/AlAlOx/Nb the correlation 
is observed in the whole temperature range. For Bi-2212 
the correlation is holding well until ^ Tc/2. The devia- 
tions at lower temperatures are caused by two artifacts 
affecting the sum-gap peak height. For Bi(Pb)-2212 the 
peaks becomes sharper because of self-heating and for 
Bi(Y)-2212 it becomes broader because of minor inho- 
mogeneity of junctions in the mesa |13j, |66j. 



B. Analysis of tunneling magnetoresistance 

Figs. [3] (a) and (b) show calculated and measured / — 
Vs of Nb/AlAlOx/Nb at T ~ 2 K for different out-of- 
plane magnetic fields. Again, a good agreement is seen. 
With increasing field, the sub-gap current increases and 
the sum-gap kink is rapidly smeared out. The I — V is 
approaching the ohmic normal-state at H ^ Hc2 ■ 

From comparison of theoretical curves in Figs. [5] (a) 
and 13] (a) it is seen that although both temperature and 
magnetic field suppress superconductivity when T — > Tc 
and H Hc2, there is a difference in how they do that. 
Namely, the temperature reduces A but does not affect 
the shape of the QP DOS, which remains gapped N{E < 
A) = and maintains a sharp BCS singularity at the gap. 
Because of that the sum-gap kink remains sharp even at 
elevated T. On the other hand, magnetic field first of all 
smears the gap singularity in the DOS and the sum-gap 



kink in I — V and increases the sub-gap conductance by 
making the DOS gapless. 

Figure [3] (c) shows simulated I — V at different H and 
at T = 0.5 K in the semi-logarithmic scale. It is seen 
that in this scale the curves remain almost parallel and 
move to lower voltage with increasing H (negative MR) . 

Fig. [3] (d) represents the similar semi-logarithmic plot 
oi I — V curves at different H and at T 2K for a larger 
Bi(Pb)-2212 mesa on the same crystal as in Fig. [5] (c). 
Here we changed the c-axis component of the field H± 
from to 17 T by rotating the crystal with respect to 
the fixed magnetic field of 17 T. Because of the extreme 
anisotropy of Bi-2212, in-plane magnetic field of 17 T 
does not produce any effect on QP DOS. This is clearly 
seen from the two rightmost I ~V, which were measured 
at H = (magenta) and at H = 17 T strictly parallel to 
the afc-planes (black). The only minor difference between 
those curves is caused by appearance of phonon-polariton 
resonances ^3], seen as small steps in I — V . 

Clearly, the general trend of experimental I — V iov 
Bi-2212 is the same as for numerical simulations and 
Nb/AlAlOx/Nb: in the semi-logarithmic scale the curves 
remain parallel and move to lower voltages, as the con- 
sequence of suppression of the superconducting gap by 
field. A similar trend was also observed for Bi-2212 in- 
trinsic Josephson junctions as a function of T [l^ [s^ j. 
For comparison, in Fig. |3] (d) we also show the I — V 
at T = 90 K~ Tc at H = 0. Apparently, the maximum 
available field of 17 T is insufficient for complete suppres- 
sion of superconductivity at T = 2 K. The parallel shift of 
I-V curves implies that the MR, V{I, H = G)-V{I, H), 
measured at fixed current and T, is approximately bias- 
independent. This is important because it allows a confi- 
dent estimation of the MR in a broad sub-gap bias range, 
including low bias with low dissipation powers, which 
precludes distortion by self- heating. 

Fig. |4] shows tunneling conductance dI/dV{V) curves 
in the semi-logarithmic scale at low T ~ 2 K and at 
different out-of-plane magnetic fields. Panels (a) and (b) 
correspond to theoretical and experimental data for the 
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FIG. 4. (Color online). Magnetic field evolution of dl /dV{V) characteristics at low T: (a) BCS calculations at T = 2 K, (b) 
the Nb/AlAlOx/Nb junction at T = 2 K and (c) a small Bi(Pb)-2212 mesa (normafized by dI/dV{T > Tc)) at T = 1.7 K. AU 
plots are shown in the semi-logarithmic scale. Low bias structure in (c) is caused by excitation of phonon-polariton resonances 
by the ac- Josephson effect [gJ . Due to state conservation, the magnetoresistance is negative in the sub-gap region, but positive 
at and above the sum-gap peak. Note also an almost parallel translation of \ndl /dV{V) curves with H. 



Nb junction from Figs. [3] (a) and (b), respectively. At 
zero field the main feature of dI/dV(V) curves is the 
sharp sum-gap peak at = 2A/e, which reflects the 
sharp BCS singularity in the QP DOS at £: = A. With 
increasing field, the peak is rapidly smeared out. Already 
at H = 0.1 Hc2 the height of the peak is suppressed 
by an order of magnitude. Simultaneously the sub-gap 
conductance at y < is growing with field. The excess 
QP current is flowing in gapless vortex cores, see Fig. [T] 
(a). It scales with the relative core area in the unit cell 
and, therefore, increases approximately linearly towards 
the normal conductance at iJ — >■ Hc2- 

In Fig. [l](c) we show dI/dV{V) characteristics in the 
semi-logarithmic scale for different H± at T = 1.7 K, for 
another mesa 1 x 1.9 /xm^ on the same Bi(Pb)-2212 single 
crystal. To simplify the analysis, the curves are normal- 
ized by the normal state curve at T = 100 K, so that the 
normal state is simply represented by the dashed horizon- 
tal line. It is seen that dI/dV{V) curves shift as a whole 
to lower voltages with increasing field. Simultaneously 
the peak looses the height while the sub-gap conductance 
vice-versa increases with increasing field. Such the be- 
havior is almost identical to that for the Nb/AlAlOx/Nb 
junction and calculations, shown in Figs. |l](a) and (b). 



gration of dl/dV curves: 

The state conservation is the reason for T- and H- 
independence of the large bias tunneling resistance i?„, 
as seen from Figs. [2] (a) and [3] (a). 

From Figs. [2] (b-d) it is clear that both in low-Tc 
and high-Tc junctions the Ohmic tunneling resistance 
Rn above the sum-gap kink is remaining almost T- 
independent. According to Eq. ([8]), this automatically 
implies that closing of the superconducting gap by tem- 
perature occurs in the state-conserving matter. 

State conservation implies that the tunneling resis- 
tance is enhanced in the sub-gap region in the same man- 
ner as the conductance is enhanced at the sum-gap peak. 
This is demonstrated explicitly, in Figs. [5] (g) and (f ) for 
Nb/AlAlOx/Nb and Bi-2212 junctions, respectively. It 
is seen that the sub-gap resistance dV/dl in the middle 
of the sub-gap region eV = A(r) and the sum-gap con- 
ductance peak grow in a similar manner with decreasing 
T . Such a scaling is an instructive way for examining the 
state conservation in Bi-2212 because it does not depend 
on a small _R„ (T) dependence 0, [s^ ' which is likely due 
to a minor in-plane (coherent) contribution to the inter- 
layer transport ^68l |. 



C. State conservation 

According to theory, suppression of superconductivity 
both by temperature and magnetic field occurs in the 
state-conserving manner. Even though superconductor- 
insulator-superconductor tunneling is not probing explic- 
itly the single QP DOS, but rather the convolution of two 
DOS from the two electrodes, Eq. ([7]), it is still possible 
to judge about the state conservation by a simple inte- 



D. The sign change of tunneling magnetoresistance 

The tunneling MR in superconducting tunnel junctions 
and the c-axis MR in layered cuprates is often described 
as negative. That is, the resistance decreases with in- 
creasing field. This is the consequence of appearance of 
the gapless state in magnetic field. The corresponding 
increase of the sub-gap DOS leads to approximately lin- 
ear increase of the low bias tunneling conductance with 
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field This is clearly seen in Fig. |4l However it is also 
seen that at the sum-gap peak the situation is reversed: 
the tunneling conductance is decreasing with increasing 
field, i.e., the tunneling MR at ^ > 1^ is positive. Thus, 
the tunneling MR is changing sign from the negative at 
the sub-gap voltage to the positive close and above the 
sum-gap voltage. Again, the behavior is similar for both 
Nb/AlAlOx/Nb and Bi(Pb)-2212 junctions. 

The sign change of the tunneling MR is a direct con- 
sequence of state conservation. The missing area of the 
sub-gap conductance with respect to normal conductance 
l/Rn is equal to the excess area of the sum-gap peak. 
Therefore, the positive MR at large bias is directly con- 
nected with the negative MR at low bias. The discussed 
sign-change of the MR is very characteristic and can be 
used for unambiguous discrimination of the supercon- 
ducting gap from non-pairing effects in the tunneling 
DOS, such as peculiarities of one-QP band structure, or 
thermal activation enhancement of the tunneling matrix 
elements for interlayer hopping [s^ j. 

V. SCALING LAWS OF TUNNELING 
MAGNETORESISTANCE 

The main purpose of this work is to determine how 
to extract useful information from tunneling MR. To un- 
derstand this we first consider the behavior of magne- 
totunneling for conventional BCS superconductors. For 
this we analyze numerically calculated and experimental 
characteristics for Nb junctions at different T and H. 



A. Zero-bias magnetoresistance 

Transport measurements are usually performed by ap- 
plying a small ac-current, i.e., probe zero-bias MR. Thick 
lines in Fig. [5] (a) show theoretical values of the zero bias 
conductance dl / dV (0) , normalized by as a function 
of H/Hc2 for T = 2, 4.7 and 7.3 K. It is seen that at low 
T = 2 K, the zero bias conductance is gradually increas- 
ing from almost zero at iJ = 0, to normal conductance 
atFc2(T) (see Fig. 11(a)). 

Temperature dependence of dI/dV{0) at H = can 
be seen from I — V curves in Fig. [2] (b). The increase of 
T leads to the decrease of A. Both factors result in the 
increased number of excited QPs above the gap, which 
initially leads to a rapid filling-in of the zero-bias dip in 
conductance and then to development of a maximum at 
V = 0. The latter represents a zero-bias logarithmic sin- 
gularity It occurs because at elevated T there is a 
substantial amount of thermally excited QPs just above 
the gap. At V = the partly filled gap singularities in 
the two electrodes are co-aligned, causing a large current 
flow from one electrode to another, which is exactly com- 
pensated by the counterfiow from the second electrode. 
However, exact cancellation is lifted at an arbitrary small 
voltage across the junction, leading to a sharp maximum 



in dl/dV. The zero-bias logarithmic singularity leads to 
an overshooting of dI/dV{0) a.t H — over the normal 
conductance at T = 7.3 K in Fig. [5] (a). 

The zero-bias logarithmic singularity makes the behav- 
ior of the zero-bias conductance non-trivial. In general, 
there is no scaling with H/Hc2 for dl ldV{Q) at different 
T, as seen from Fig. [S](a). To avoid complications caused 
by the zero-bias singularity we look at the behavior of the 
conductance at finite bias. 



B. Scaling of the sub-gap conductance and the 
sum-gap resistance 

Fig. [5] (b) shows field dependence of the sub-gap con- 
ductance at the middle point eV = A for the same tem- 
peratures as in (a). It is seen that unlike dl /dV{Q), the 
sub-gap conductance is showing a fairly universal linear 
scaling as a function of H/Hc2 for different T. 

As follows from Figs. [2](g) and (h), state conservation 
implies that the deficit of the sub-gap conductance is 
directly connected to the excess of the sum-gap peak, i.e., 
the deficit of the sum-gap resistance. Fig. [5] (c) shows 
magnetic field dependence of the sum-gap resistance. It 
is showing an almost universal, slightly non-linear scaling 
as a function of H / Hc2 in the wide T-range. 

Fig. [5](d) shows magnetic field dependence of the sum- 
gap peak voltage normalized by that at zero field. It 
exhibits a very simple universal linear scaling as a func- 
tion of H/Hc2 in the whole T— range. Interestingly, the 
peak voltage does not go to zero at H = 11^2 , but rather 
stops half-way at eV ~ A{H = 0). This is due to an 
interplay between the reduction of A(i/), see Fig.[T](b), 
which moves the peak down, and simultaneous strong 
smearing of the maximum in the QP DOS, see Fig.[T](c), 
which moves the peak up in voltage. At H = Hc2, A 
does vanish, but voltages/energies of very broad maxima 
in dl/dV or spatially averaged DOS do not vanish (4^ . 
as can be seen from Fig. H (b). 



C. Extraction of Hc2 for Nb 

The universal scaling of the sum-gap peak voltage 
Vp (H/ Hc2) can be used for extraction of Hc2 from magne- 
totunneling data. Symbols in Fig. [S] (d) show the results 
of fitting for the Nb/AlAlOx/Nb junction, using Hc2{T) 
as the only adjustable parameter. It is seen that the 
agreement with theoretical calculations is excellent. T- 
dependence of the obtained upper critical field is shown 
in Fig. [5] (e). Due to the proximity effect in Al it is 
more linear at low T than that for pure Nb. A similar 
H^{T) was also reported for proximity coupled Nb/Cu 
multilayers [TOj . 

Using thus extracted Hc2, we check the scaling of other 
experimental parameters such as the zero bias and the 
sub-gap conductances and the sum-gap resistance, shown 
by symbols in Figs. [S] (a), (b) and (c). In all cases the 
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FIG. 5. (Color online). Analysis of scaling of magnetotunneling characteristics at different T as a function of H / Hc2 in theory 
(thick lines) and for the Nb junction (symbols), (a) Zero bias conductance. Note that there is no scaling of dI/dV(Q) for 
different T due to progressive development of the zero-bias logarithmic singularity with increasing T. Thin solid lines and the 
inset show evolution of the conductance of Nb-electrodes at Ji" > Hc2- It is seen that the surface superconductivity in electrodes 
survives up to Hd — 1.69//c2- (b) Scaling of the sub-gap conductance at finite bias V = A/e. (c) Scaling of the sum-gap 
resistance, (d) Scaling of the sum-gap peak voltage, (e) The upper critical field for Nb, extracted from the scaling. 



agreement with theory is very good, confirming the cor- 
rectness of determination of iJc2 • 



D. Surface superconductivity above Hc2 in Nb 

The only case where there is a certain disagreement be- 
tween theory and experiment is dl /dV{Q) at high fields 
in Fig. [5] (a): experimental data at T = 2 K and 4.7 K is 
clearly not reaching the normal conductance at H = Hc2 ■ 
At higher fields it is not possible to measure junction 
characteristics because the Nb electrodes are no longer 
capable to carry the supercurrent and turn into the resis- 
tive state. Since the electrode resistance is several times 
larger than that of the junction, the measured resistance 
in this case is mostly given by the longitudinal resistances 
of the Nb electrodes. This leads to the drastic decrease 
in the measured zero-bias conductance, as shown by thin 
solid lines in Fig. [S](a). From this data it is obvious that 
some residual superconductivity in the Nb electrodes is 
remaining up to fields significantly larger than Hc2- 

The inset in Fig. [S] (a) shows the detailed view of 
the onset of the superconducting transition in Nb elec- 
trodes. It is seen that at low T = 2 K superconductivity 



in Nb survives up to almost exactly 1.69-ffc2, which is 
the expected value of the third critical field iJc3- At 
Hc2 < H < Hc3 superconductivity exists only in surface 
layers. A similar behavior has been reported in clean Nb 
[tH - ItsI as well as in MgB2 (tJ- From the inset in Fig. 
[S] (a) it is seen that the ratio i?c3 / Hc2 is decreasing with 
increasing temperature. A similar behavior was reported 
for clean Nb and discussed in terms of a tricritical point 
(rij . However, it should be noted that our Nb films are 
in the dirty limit and are affected by the proximity effect 
with Al. The columnar structure of sputtered Nb films 
with a large effective surface-to- volume ratio and columns 
orientation perpendicular to the film may also enhance 
the role of surface superconductivity in the out-of-plane 
magnetic field. 



VI. ANALYSIS OF INTRINSIC 
MAGNETOTUNNELING IN BI-2212 

Observation of a fairly conventional behavior of tun- 
neling MR in overdoped Bi(Pb)-2212 encourage us to 
employ the derived scaling laws for extraction of the T- 
dependent upper critical field, which remains a contro- 
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FIG. 6. (Color online), (a) Field-dependence of the sum-gap resistance in Bi(Pb)-2212 at different T (normalized to the normal 
state value). Dotted lines indicate extrapolations to H = Hc2- (b) Magnetic field evolution of dl /dV{V) curves at T = 75 K 
(normalized by dI/dV{V) at T ~ Tc). (c) Relative change of dl /dV upon rotation of the crystal in field // = 10 T from the 
in-plane to the c-axis direction. Inset shows enhanced view of the curve at T = 95 K. 



versial issue for cuprates, as mentioned in the introduc- 
tion. The obtained scahng of the sum-gap peak voltage 
Vp{H/ Hc2) is vaUd not only for tunnel junctions made 
of s-wave superconductors. A similar scaling was also 
reported for A(H/Hc2) in the mixed state of c?-wave su- 
perconductors !51|. In the remaining part of this work we 
apply the scaling rules for Bi(Pb)-2212 mesas in order to 
understand how normal or abnormal the behavior of in- 
trinsic tunneling magnetoresistance is and in an attempt 
to estimate the upper critical field. 



A. Scaling of the sum-gap peak 

As seen from Fig. S] (c), at low T dI/dV{V) char- 
acteristics of the Bi(Pb)-2212 mesa, normalized by the 
normal state dl /dV , behaves in a conventional manner: 
the sum-gap peak is smeared and moves to lower voltage 
with increasing the c— axis component of the field. Solid 
squares in Fig. [5] (a) represent the corresponding peak 
resistance as a function of field. Apparently, it follows the 
same tendency of quasi-linear grows as for Nb-junction, 
shown in Fig. [S] (c). This allows an approximate estima- 
tion of Hc2 by linear extrapolation of the curves to unity, 
as shown by dotted lines. 

Figure H] (b) shows magnetic field dependence of 
dI/dV{V) for the mesa from Fig. |3](d) at T = 75 K. Here 
we also normalized the dl/dViV) curves by that above 
Tc- It is seen that the general behavior of the tunneling 
MR is the same as at low T: the sub-gap conductance 
increases (negative MR) at the expense of the sum-gap 
conductance peak (positive MR) in the state-conserving 
manner. The corresponding peak resistances are shown 
by open circles in Fig. in](a). However, the sum-gap peak 
is no longer moving to lower voltages with increasing H , 
but instead spreads out to higher voltages. As discussed 



in Ref. [4J], such a behavior is not totaly unusual. As a 
matter of fact the upturn of the peak voltage for high H 
and T is also observed in theory and for Nb/AlAlOx/Nb 
junctions, even though in a smaller scale, see curves at 
T = 4.7 K in Fig. [S] (d). As discussed above, the out- 
ward motion of the peak is the consequence of smearing 
of the peak in DOS, rather than the actual increase of 
A. This is also the reason why Vp does not go to zero at 
H Hc2, despite A does vanish. 

Figure [5] (c) shows the ratio of conductances aX H = 
10 T oriented parallel to a6-planes (which, as demon- 
strated in Fig. [21(d) is equivalent to zero- field) and in the 
c-axis direction. Such a normalization perfectly removes 
all field independent features and allows observation with 
unprecedented clarity superconducting parts of the spec- 
tra at r — i> Tc- In Fig. [5] (c) we clearly see the sum-gap 
peak at T = 85 K, which is practically indistinguishable 
in dIldV{V) characteristics. Remarkably, we can ob- 
serve the peak even at T = 95 K (see the inset), which is 
above the critical temperature for appearance of phase- 
coherence in the c-axis direction J'p''"*'^ = 91 _ 92 K. 
Apparently the superconducting gap is still present at 95 
K, but its value is rapidly decreasing at this temperature, 
as it does in BCS theory close to the mean- field Tco, see 
Fig. [5] (f). This is consistent with the conclusion of Rcf. 
[13i] that superconductivity in near optimally doped Bi- 
2212 appears by means of the second-order phase transi- 
tion in the conventional BCS-manner. However, the true 
thermodynamic mean-field critical temperature Tco — 96 
K is some- what higher than T^'^"'^'^. 



B. Scaling of the sub-gap voltage 

As seen from Fig. [3](c), the voltage at a given current 
is decreasing with increasing field and reaches the normal 
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FIG. 7. (Color online). Magnetic field dependence of the voltage at a given current for (a) theoretical calculation, current 
levels correspond to horizontal lines in Fig. |3](c), and for Bi(Pb)-2212 mesa at (b) T = 3 K and (c) T = 70 K at different bias 
currents. Inset in (b) shows normal-state / — Vs obtained by adjustment of V — Ki so that they vanish at the same H = Hc2, 
as indicated by dotted lines in (b) and (c). 



state value Vn{I) = IRn at H = Hc2- Fig. [7] (a) shows the 
corresponding values of V{I,H) for four bias currents, 
indicated by horizontal lines in Fig. [3] (c) , as a function 
of H/Hc2 ■ It is seen that in a wide bias range (an order 
of magnitude) V{I,H) — IRn have the same magnitude 
and disappears in a linear manner at — > Hc2- This 
fairly good and almost bias-independent scaling is the 
consequence of the almost parallel shift of / — T^s (in the 
semi-logarithmic scale) with increasing as shown in 
Figs. |3](c), . The dI/dV{V) curves behave in the same 
manner, as shown in Fig. 21 

Fig.|3](d) demonstrates that for Bi(Pb)-2212 such par- 
allel translation of I — V curves is even more impressive 
and expands to almost two orders of magnitude in bias 
current. However, there is one major obstacle for deter- 
mination of the Hc2 from such the scaling. Namely, the 
shape of the / — in the normal state, i.e., with com- 
pletely suppressed superconductivity but at low T < Tc, 
is unknown. We can only say for sure that it remains 
non-linear above Tc as a result of thermal-activation c- 
axis transport (not necessarily connected with the pseu- 
dogap) [34]. To go around this problem, we note that 
according to Fig. [7] (a), voltages at all bias levels must 
reach the normal state values Vn{I) in a quasi- linear man- 
ner at the same H = Hc2- Therefore, we used Vn{I) as 
an adjustable parameter for each bias current, so that 
all V{I) — Vn{I) voltages go to zero at the same point 
H — Hc2{T), as indicated by dotted lines in Figs. [7] (b) 
and (c). Thus obtained normal state I — V curves at 
T < Tc are shown in the inset of Fig. [7] (b). 

Figure [7] (b) and (c) show i?j_-dependence of the sub- 
gap voltages at different bias currents, for (b) T = 3 
K and (c) 70 K. We observe a quasi-linear reduction of 
voltages, consistent with theoretical curves in Fig. [7] (a). 
The extrapolated value of Hc2 at T = 3 K is 75.5 ± 
14. 5T. A significant uncertainty is due to the need for 



remote extrapolation from the maximum available field 
of 17 T. However, at low T it is not crucially affected by 
the adjustment of Vn- Indeed, already from the raw data 
in Fig. [3] (d) it is seen that in fields from to 17 T the 
I—V curves made approximately 1/4 — 1/5 of the journey 
to the normal state I —V at T = 90K. Assuming a linear 
dependence of V{H), this provides a similar estimation 
of Hc2 ~ 68 - 85 T. 



C. Temperature dependence of the maximum 
c— axis magnetoresistance 

Fig. [5] (a) shows temperature dependence of the max- 
imum measured shift of voltage at a constant current 
upon variation of the c-axis magnetic field from to 17 
T. It is seen that the negative sub-gap tunneling MR 
is rapidly decreasing with increasing temperature, cf. 
scales in Figs. [7] (b) and (c). This makes determina- 
tion of Hc2 at elevated T less confident because it be- 
comes more sensitive to exact values of Vn ■ For compar- 
ison, Fig. [5] (b) shows similar data for the Nb-junction 
at H ^ 0.44 T, which corresponds to the same ratio 
H/Hc2{T = 0) ~ 0.24 as for Bi(Pb)-2212 from panel 
(a). The tunneling MR in the Nb-junction vanishes at 
T ~ 7 K< Tc, at which i/c2(T) = H, see Fig. [5] (e). 

From comparison of Figs. [8] (a) and (b) it is seen that 
behavior of the maximum MR is qualitatively similar for 
Nb and Bi(Pb)-2212. However, for Bi(Pb)-2212 it be- 
comes more smeared and fuzzy at T Tc. To some 
extent the smoother T— dependence for Bi-2212 may be 
caused by the d-wave symmetry, because nodal QP's are 
more prone to thermal activation above the gap at ele- 
vated T. Still this does not explain everything. 

The inset in Fig. [5] (a) shows the same data for Bi(Pb)- 
2212 in the semi-logarithmic scale. It is seen that the sub- 



12 




O 
< 



0,0 



40 80 120 

Bi(Pb)-2212 



1 I 1 Wa-ftn ^ rin 



40 80 

r(K) 



120 



0,8 - 



> 



0,6 - 



B 

^ 0,4 

I 

O 



0,2 - 



0,0 



-I — I — I — I — I — I — I — I — J- 

... Nb/AlAlOx/Nb 



-(b) 



H = 0.44 T 

_i I 1 I 1 L 



4 6 

r(K) 



80 



P 60^ 1 

ti::"4o- (c) 



20 - 



r'''''''=91K Bi(Pb)- 



I 



From: 



sub-gap V{H) 
O sum-gap dV/d/ 



2212 - 




i r =96K. 

cO 



40 80 

T(K) 



120 



FIG. 8. (Color online). Temperature dependence of the maximum magnetoresistance at a constant current, AV(0 — H) — 
V{I,H = 0) - V{I,H), (a) for Bi(Pb)-2212 mesas at // = 17 T, (b) for the Nb/AlAlOx/Nb junction at ~ 0.44 T with the 
same H/Hc2{0) ratio. The inset in (a) shows the same data in the semi-logarithmic scale. It is seen that the negative sub-gap 
magnetoresistance decays almost exponentially with increasing T, but retains a finite value at T > Tc. (c) Extracted upper 
critical field in the c-axis direction for Bi(Pb)-2212. 



gap MR start to rapidly decrease (approximately expo- 
nentially) with increasing temperature at 40 K< T < Tc 
(as indicated by the dotted line), experiences an addi- 
tional drop at Tc (marked by the dashed vertical line) 
and then continued to decrease at even faster rate above 
Tc- We assume that the remaining small negative MR 
above Tc is a consequence of fluctuation superconductiv- 
ity, which becomes progressively less sensitive to mag- 
netic field due to growth of the effective H*2 (x T — Tc 
above Tc |29|]. 



D. Extraction of Hc2 for Bi-2212 

As follows from Fig. [5] (c), scaling of sum-gap magne- 
toresistance provides the most accurate way of determi- 
nation of Hc2- The corresponding data for Bi(Pb)-2212 
are shown in Fig. [n](a). Data for 1.7 K and 75 K corre- 
spond to Figs. |4](d) and[6](b). The upper critical field is 
estimated using a linear extrapolation towards the nor- 
mal resistance, as indicated by dotted lines. There is a 
significant uncertainty in such extrapolation, however it 
is more robust than that made from sub-gap MR, see 
Fig. [71 because the normal resistance at a large bias is 
more unambiguous and has a weak T-dependence 0, [3^ , 
as seen from Fig. [5] (c). 

Figure [5] (c) shows the extracted T-dependence of the 
upper critical field for the slightly overdoped Bi(Pb)-2212 
single crystal. Solid squares and open symbols are ob- 
tained from analysis of scaling of the sub-gap voltage. 
Figs. [3 (b) and (c), and the sum-gap resistance. Fig. [S] 
(a), respectively. Hc2{T = 0) is ~ 70 T and certainly 
decreases with increasing T. As mentioned in the intro- 
duction, one of the principle questions is whether Hc2 
goes to zero at Tc- This appears to be a difficult ques- 
tion. First, the notion of Tc is fuzzy, as indicated by the 



gray area in Fig. [51 We see the clear presence of the su- 
perconducting gap in the DOS above the phase coherent 
rj^phase ^ 91 K up to the mean-field T^o ~ 96 K jl3|. The 
main experimental challenge is associated with a very 
small sub-gap MR close to T^ see Fig. [SI (a). However, 
at T = 95 K any sign of the sum-gap peak is absent at 
H = 15 T, which provides a rough estimation for the last 
point in the Hc2{T) diagram of Fig. [8] (c). 



VII. CONCLUSIONS 

We performed a detailed comparative analysis of 
tunneling magnetoresistances in conventional low-Tc 
Nb/AlAlOx/Nb junctions, small Bi(Pb)-2212 intrinsic 
Josephson junctions and microscopic calculations. It 
was found that magnetotunneling in slightly overdoped 
Bi(Pb)-2212 is qualitatively similar to that in conven- 
tional BCS-type superconductors. 

From the data presented above it is clearly seen that 
both temperature and magnetic field suppress supercon- 
ductivity in Bi(Pb)-2212 in the state conserving manner. 
Magnetotunneling provides a particularly clear demon- 
stration of this: due to conservation of states the MR 
changes sign from the negative in the sub-gap, to the 
positive at the sum-gap bias. Continuing strong MR 
well above the sum-gap peak with powers up to several 
times that on the peak indicate that the mesa remains 
in the superconducting state and the peak is not caused 
by self-heating. This clearly demonstrates that intrinsic 
tunneling can provide unambiguous information about 
bulk electronic spectra of Bi-2212. 

Observation of state conservation implies that QP 
states released upon suppression of superconductivity by 
magnetic field are not taken over by a competing order, 
like charge or spin density wave. In other words, there 
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is no field-induced non-superconducting order in slightly 
overdoped Bi(Pb)-2212. The situation may, however be 
different for underdoped Bi-2212, for which non-state- 
conserving characteristics have been reported [1^, [t^ . 

We derived theoretically and verified experimentally 
scaling laws of various magnetotunneling parameters. 
Those scaling laws were employed for accurate extraction 
of the upper critical fields and in the case of Nb provided 
a clear evidence for the existence of an extended region 
of surface superconductivity at Hc2 < H < Hcs- 

For Bi(Pb)-2212, it was found that Hc2{T = 0) ~ 
70 T and decreases significantly upon approaching Tc- 
The parameters of Bi(Pb)-2212 were obtained from self- 
consistent analysis of magnetotunneling data at different 
levels of bias, dissipation powers and different mesa sizes, 
which precludes the influence of self-heating. The ampli- 
tude of the sub-gap magnetoresistance is suppressed ex- 



ponentially at T > Tc/2. It remains negative, although 
very small, above T^, probably indicating existence of an 
extended fluctuation region. 

We conclude in general that intrinsic magnetotunnel- 
ing in small mesa structures is a very powerful tool for 
analysis of bulk superconducting properties of cuprates. 
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